function crit=estimMonteCarlo20191030(theta)
global Glob
dimNation=Glob.dimNation;
dimProv=Glob.dimProv;
TrueMarriage=Glob.TrueMarriage;
omega=theta(1:3);
sigma=theta(4);
alpha=repmat(eye(dimNation),[1 1 dimProv]);
ind=logical(1-eye(dimNation));
ind=repmat(ind,[1 1 dimProv]);
alpha(ind)=theta(5:end);

param.alpha=alpha;
param.omega=omega;
param.sigma=sigma;

PredMarriage=solveMonteCarlo20191030(param,Glob);

Diff=PredMarriage-TrueMarriage;
Diff=Diff(:);
crit=Diff'*Diff;